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We given an exposition of a proof that a mirror symmetric config- 
'— uration of two convex analytic obstacles in (M) 2 is determined by its 
. Dirichlet resonance poles. It is the analogue for exterior domains of 

| the proof that a mirror symmetric bounded simply connected analytic 

plane domain is determined by its Dirichlet eigenvalues. The proof 
. uses 'interior/exterior duality' to simplify the argument. 

(N : 
o . 

§ : 1 Introduction 

This article is part of a developing series [Zl, Z2] concerned with the inverse 
spectral problems for analytic plane domains. It is essentially the lecture 
we presented at the IMA workshop on Inverse Spectral Problems in July, 
| 2001 and represents the state of our knowledge, methods and results at that 

time. We apply these methods to prove an analogue for exterior domains of 
the result proved in [Z2] for interior ones, namely that a mirror symmetric 
configuration of two convex analytic obstacles in (R) is determined by its 
Dirichlet or Neumann resonance poles. The proof combines the known result 
that wave invariants of an exterior domain are resonance invariants (see §3) 
with the method of [Zl, Z2] for calculating the wave invariants explicitly 
in terms of the boundary defining function. In keeping with the expository 
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nature of the lecture, we give a detailed exposition of the background results 
on the Poisson relation for exterior domains and of the main steps in [Zl, Z2] 
in the calculation of wave invariants (§6) and on determining the domain 
from its wave invariants (§7). 

The motivating problem is whether analytic domains are determined by 
their spectra for Dirichlet or Neumann boundary conditions. This inverse 
spectral problem can be posed for both interior and exterior domains. For 
the interior problem, we assume the domain is a bounded simply connected 
plane domain and ask whether its Dirichlet or Neumann spectrum determines 
the boundary up to rigid motion. For the exterior problem, we assume the 
domain is the complement of two bounded, simply connected obstacles and 
ask whether we can determine the pair of obstacles from the resonance poles 
(or scattering phase) of the exterior Dirichlet or Neumann Laplacian. Since 
the method is similar for both boundary conditions, we assume the boundary 
conditions are Dirichlet in the exterior domain. 

In the previous paper [Z2], the interior inverse spectral problem is studied 
for analytic domains with a certain mirror symmetry a. The symmetry is 
assumed to fix a bouncing ball orbit (as a set) and to reverse its orientation; 
i.e. if the domain is translated and rotated so that the bouncing ball orbit lies 
along the y-axis with its midpoint at the origin, then a(x,y) = (x, —y). In 
addition, we assume the length of the bouncing ball orbit is a fixed number L 
and that the orbit satisfies a non-degeneracy condition. The main result (see 
Theorem (1.4)) is that the interior domain satisfying some generic conditions 
is determined by its Dirichlet spectrum among other such domains. We note 
that this result is stronger than the one in [Z3] in having eliminated one of the 
two symmetries. We should emphasize that the proof does not immediately 
imply the result for the other mirror symmetry a(x, y) = (—x, y), which fixes 
7 pointwise and preserves the orientation. 

Our main purpose in this paper is to extend the result to connected 
exterior domains of the form 

n c ■= (R) 2 - n, with n = {o u t x , l O) (i) 

where O is a bounded, simply connected analytic domain and where 

t~x,l = reflection through the line (i/ x , x + l/2Lv x — y) = 0. 

Here, x e O, v x is the outward unit normal to O at x and L > is a 
given positive number. Thus, the obstacle consists of two (non-intersecting) 
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isometric components which are mirror images of each other across a common 
orthogonal segment of length L. Such an obstacle is called a (symmetric) 2- 
component scatterer. The segment is the projection to (M) 2 of a bouncing 
ball orbit 7 of the exterior billiard problem. Throughout, the notation Q 
refers to a bounded domain and Q c is the unbounded complement. 



The problem is thus to recover the obstacle O U t x ^O from the set of 
resonance poles of the exterior Dirichlet Laplacian A^> c , i.e. the poles 
of the analytic continuation of its resolvent 



to the logarithmic plane. To be precise, we consider the class OBSTA- 
CLES^ of two-component obstacles satisfying: 




Figure 1: Mirror symmetric 2-obstacle scatterer. 
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• (i) O is simply connected and real analytic; 

• (ii) 7 is a non-degenerate bouncing ball orbit, whose length L is isolated 
in LspiVl) U Lsp{Vt c ). 

We denote by TZES{Vt c ) the set of resonance poles of the Laplacian A^, 
of the domain Q c with Dirichlet boundary conditions. 
We now state the main results. 

Theorem 1.1. TIES : OBSTACLE^ h-> C N is 1-1. 

In fact, in combination with a result of M. Zworski [Zw], the proof shows 
more: the obstacle is determined by the resonances close to the real axis 
which are associated to the bouncing ball orbit 7 of the exterior billiard 
problem. We note that assumption (ii) is stronger than the one in the interior 
case (cf. Theorem (1.4) in demanding multiplicity one for L in the combined 
interior and exterior length spectra. By using the original proof in [Zl, Z2], 
one could remove Lsp(Q) from the assumption and just demand multiplicity 
one in Lsp(Q c ). We do not do so, because we wish to present a simpler proof 
than the one in [Zl, Z2], as we will explain below in the introduction. 



As a corollary, we prove that the obstacle is determined by its scattering 
phase. We denote by 5d(A) the Dirichlet scattering operator for Q c and by 
«d(A) = det Sd(A) the scattering phase. 

Corollary 1.2. // two exterior domains Q in OBSTACLEi 2)L have the 
same scattering phase sd(A), then they are isometric. 

We note that the length L of the segment between the components has 
been marked, since apriori it is not a resonance invariant. Under some ad- 
ditional assumptions, it is a resonance invariant and one does not have to 





Figure 2: Bouncing ball orbit. 
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mark it. For instance, if one assumes that O is convex, then the bouncing 
ball orbit 7 is the unique periodic reflecting ray of the exterior domain up 
to iterates, so its length is a resonance invariant. Thus, let COBSTACLE% 2 
denote the convex component obstacles in OBSTACLE^. We then have: 

Corollary 1.3. TIES : COB STAC LE% 2 i-> C N is 1-1. 

For the sake of completeness, we also recall the precise statement of the 
interior result of [Z2]. Let Spec(f2) denote the spectrum of the interior Lapla- 
cian of the domain Q with either Dirichlet or Neumann boundary condi- 
tions. 

Theorem 1.4. Let T>l denote the class of simply connected real-analytic 
plane domains Q satisfying: 

• (i) There exists an isometric involution a ofQ; 

• (ii) a 'reverses' a non- degenerate bouncing ball orbit 7 — > 7 _1 ; 

• (Hi) The lengths 2rL of all iterates Y have multiplicity one in Lsp(Q), 
and the eigenvalues of the linear Poincare map P 7 are not roots of 
unity; 

Then: Spec: T> 1L \— > (M) 1 ^ is 1-1. 

To be cautious, we should point out that these results are part of a work- 
in-progress which has not yet reached its final form. They have been inde- 
pendently verified when there is an additional mirror symmetry (i.e. each 
obstacle is itself 'left-right' mirror symmetric), and in this case the relevant 
calculations of wave invariants are easy. When there is only one mirror sym- 
metry, the calculations become somewhat messy (see [Z2]) and the proof 
therefore becomes rather 'unstable'; at this time of writing, some details 
have been independently verified (by R. Bacher [B]), but some have not. 

Let us briefly describe the proof and the organization of the paper. As 
emphasized by Zworski, there is nothing really new in Theorem 1.1 beyond 
the result of [Zl, Z2] in the interior case (i.e. Theorem 1.4). We are just 
combining the known fact that wave trace invariants of the exterior domain 
are resonance invariants with the calculation of the wave trace invariants in 
terms of the germ of the defining function of dQ at the endpoints of the 
bouncing ball orbit. This calculation was done for the germ of any kind of 
bouncing ball orbit in [Z2], developing a method originating in the work of 
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Balian-Bloch [BB1, BB2]. However, it seems to us worthwhile to collect all 
the facts one needs for the proof of Theorem (1.1) in one place and to explain 
the main steps without all the details in [Zl, Z2]. 

Moreover, our exposition has one novel point: in §4, we use the so-called 
interior /exterior duality to simplify formula for the resolvent trace in [Zl]. 
Interior/exterior duality is a distributional trace formula, which we informally 
write for r >> as 

Tr \Rf[R% C (k + it) © R%(k + it) - R (k + ir)] 

d ( \ (2) 

= — logdet (l + N(k + iT)^ 
dk \ 

Here, the determinant is the Fredholm determinant, R^, Rp denote Neu- 
mann (resp. Dirichlet) resolvents on X and the notation Trx indicates the 
space on which the trace is taken. We write L 2 ((R) 2 ) = L 2 (Q) © L 2 (Q C ) and 
view R^(k + it) © R%(k + it) as an operator on this space. For notational 
simplicity, we do not put in the explicit projections 1^, Iqc (i.e. the charac- 
teristic functions). Also, N(k + it) is a boundary integral operator which 
will be defined in §4. The correct statement and proof of (2) will be given in 
Proposition (4.1). It is essentially a differentiated version of the formula 

^ log det S D (E) + ^N N (E) = ^ log Cd,n(E + <0) (3) 

of Eckmann-Pillet (cf. [EP, EP2]), relating the scattering phase of the exte- 
rior Dirichlet problem, the eigenvalue counting function of the interior Neu- 
mann problem and a spectral zeta function for the integral operator iV along 
the boundary. 

This reduction to the boundary trace on the right side of (2) simplifies 
some of the technical details of [Zl] and allows one to deduce the inverse 
spectral result simultaneously for the interior and exterior problems. In ad- 
dition, it brings our calculations into closer contact with the physics litera- 
ture, where the Balian-Bloch approach is now almost always applied to the 
boundary trace (see e.g. [AG, THS, THS2, GP]. The price we pay is that 
the combination of inside and outside requires us to make the additional 
multiplicity free assumption on the combined inner and outer length spectra 
to obtain the inverse results in Theorem 1.1. As mentioned above, this as- 
sumption could be eliminated; but we feel that the simplification in the proof 
is worth the extra assumption. To simplify the exposition, we also assume 
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in the last step that O is convex. Hence, we only complete the proof for 
hyperbolic bouncing ball orbits. For the general case, we refer to [Z2]. 

Let us give a brief outline of the method of [Zl, Z2] and this paper for 
determining an analytic domain from its spectrum. It is based on the use of 
certain spectral (i.e. resonance) invariants known as the wave trace invariants 
at the bouncing ball orbit. As the name implies, the wave trace invariants 
at a periodic reflecting ray 7 are coefficients of the singularity expansion of 
the relative trace of the wave group E^°(t) = costA/A^ c at t = L 7 (the 
length). The inverse results are proved by explicitly calculating the wave 
trace invariants at a bouncing ball orbit in terms of the Taylor coefficients of 
the boundary defining function at the endpoints, and then determining the 
boundary from these Taylor coefficients. 

Rather than studying the wave trace per se in [Zl, Z2] we follow [BB1, 
BB2] in studying an essentially equivalent expansion involving the resolvent. 
Potential theory gives an exact formula for the Dirichlet (or Neumann) re- 
solvent of a domain in terms of the free resolvent of (R) 2 and of the layer 
potentials associated to the domain. Formally, one may derive an infinite 
series expansion known as the 'multiple reflection expansion' of [BB1, BB2], 
whose Mth term corresponds intuitively to M reflections on the boundary. 
The trace of each term is, again formally, an oscillatory integral correspond- 
ing to M-fold periodic reflecting rays. 

The advantage of this approach is that the formula is exact and does not 
require a microlocal paramatrix construction, which is messy and complicated 
for bounded domains and therefore very hard to use in inverse spectral theory. 
Moreover, the terms of the mutiple reflection expansion are 'canonical' in 
that the integrands are the same for all domains; the only difference lies in 
the domain one integrates over. Once it is 'legalized', it gives an explicit 
and even routine algorithm for calculating wave invariants. At this time 
of writing, no other approach seems to provide a workable algorithm for 
doing the calculations. The disadvantages are that the legalization has two 
complicated parts: one needs to estimate the remainder in the infinite series 
and to regularize the oscillatory integral defined by the M term. 

Our formulae for the wave invariants come from applying stationary phase 
to the regularized traces. One of the principal results is that, for each iterate 
7 r , only the term M = rm of the multiple reflection expansion contributes 
important inverse spectral data, because only this term contains the max- 
imum number of derivatives of the defining function of Q at a given order 
of k~i . The heart of the matter is then the calculation and analysis of the 
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wave invariants. Since one needs wave invariants of all orders, one seemingly 
runs into an infinite jungle of complicated expressions. To tame the jun- 
gle, we enumerate and evaluate the terms of the wave trace invariants using 
Feyman diagram techniques. It turns out that only five diagrams play an 
important role (this has been confirmed by R. Bacher). The remaining issue 
is the calculation of the amplitudes and their dependence on r. It is this 
latter dependence that is crucial in allowing one to remove a symmetry of 
the domain. 

The calculation of wave trace invariants is the same for periodic reflecting 
rays of interior and exterior domains since it depends only on the germ of 
the boundary at the reflection points. The material in §5 - §7 summarizes 
the calculation in [Zl, Z2]. We hope that this guide clarifies the calculation 
and makes it easier to check the details. 

We close the introduction with some remarks on open problems in the 
inverse spectral problem for analytic plane domains. The most obvious one 
is whether one can eliminate the remaining symmetry assumption. It may be 
that one can recover the domain from wave invariants at one closed billiard 
orbit, or that one has to combine information from several orbits. A key 
problem is, how much of the Taylor expansion of the boundary defining 
function one can recover from wave invariants at one orbit? 

The method of this paper and [Z2] gives detailed but somewhat undi- 
gested information about wave invariants in terms of the Taylor coefficients 
of the boundary defining function at endpoints of a bouncing ball orbit. An- 
other approach, suggested originally by Colin de Verdiere [CdV], is to focus 
on the Birkhoff normal form of the billiard map at the closed orbit. As was 
proved by V. Guillemin in the boundaryless case [G], the Birkhoff normal 
form at a closed orbit is a spectral invariant. In [Z3] the author generalized 
this to the boundary case. Colin de Verdiere's observation that the Birkhoff 
normal form at a bouncing ball orbit determines an analytic plane domain 
with two symmetries then gave the solution of the inverse spectral problem 
for analytic plane domains with these symmetries [Z3]. We found it difficult 
however to calculate the normal form coefficients by this method without 
assuming two symmetries. That motivated us to try the Balian-Bloch ap- 
proach. 

In the recent paper [ISZ], Iantchenko-Sjostrand-Zworski give another proof 
of the inverse result with two symmetries using Birkhoff normal forms of the 
billiard map and quantum monodromy operator rather than the Laplacian; 
the method is quite elegant and flexible. However, we believe that their nor- 
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mal form of the monodromy operator is the same as the normal form in [Z3] of 
what is called there the semiclassical wave group. We changed to the present 
approach because we found it too difficult to calculate the latter without 
the two symmetry assumption. However, the possibility remains that the 
Birkhoff normal form could be a better way to organize the inverse spectral 
data than the wave invariants themselves. The question is, how much of the 
boundary defining function is determined by the Birkhoff normal form of the 
wave group or the monodromy operator at a periodic orbit? It is known that 
the classical Birkhoff normal form by itself will not determine the Taylor 
coefficients of the domain unless there are two symmetries. 

Acknowledgements. We are heavily indebted to M. Zworski for many 
helpful conversations on the contents of this article, particularly during a 
visit at the Erwin Schrodinger Institute during June 2001. In addition, we 
thank the ESI for making the visit possible. We also thank J. Ralston for 
confirming that various results of [BGR] which we need are also valid in the 
even- dimensional case. Finally, we thank the IMA for the opportunity to 
present the results in their final form. 

2 Billiards and the length functional 

Throughout this paper, we will assume that our obstacle is up-down symmet- 
ric in the sense that there is an isometric involution a of Q which interchanges 
the endpoints of a bouncing ball orbit (extremal diameter). However, we only 
use the assumption in the last step of calculating wave invariants and deter- 
mining the domain. Elsewhere, we make it for notational convenience. We 
now extend the notation and terminology of [Zl] so that it applies to exterior 
domains which are complements of obstacles with two components. 

We align the obstacles so that the bouncing ball orbit 7 projects to the 
line segment a_a + along the y-axis and so that its midpoint lies at the orgin 
G (1R) 2 and its endpoints are a± = (0, ±L/2). We refer to the top component 
of Q as the + component + and the bottom one as the — component O-. 

There are two convenient parametrizations of dfl. First, we denote by 
q±(9) the arc-length parameterization of 0± with q±(0) = a±. At the end we 
will also use graph parametrizations: In a small strip T e (j) around aia 2 , the 
boundary consists of two components which are symmetric graphs over the 
x-axis. We write the graphs in the form y = ±f(x) near a±. 

In the multiple reflection expansion of the exterior resolvent R^{k + 
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ir), we will encounter multiple integrals f^ dn ^ M - Since dfl consists of two 
components, each such term breaks up into 2 M multiple integrals over the 
circle T. The terms can be enumerated by maps a : {1, . . . , M} — > {±}. We 
needed the same enumeration in [Z2] to denote the two local components of 
Q at the ends of the bouncing ball orbit. 



2.1 Length functional 

For each map a : {1, . . . , M} — > {±} as above, we define a length functional 
on T M by: 



L CT ((/9l, . . .,(p M ) 

= kfliW - *x 2 (y?2)| + • • • + (g^.^M-i) - q* M (¥M) 



(4) 



It is clear that L CT is a smooth function iia(j) ^ a(j + l) for any j. When 
there do exist such j, then L a is singular on the 'large diagonal' A,j+i : = 
{<fj = (pj+i}, where it has singularities . A standard (and easy) calculation 
shows that 

d 

At = sin Z(q a {(p j+1 ) - q aj ((pj), )) 



- sin (</?,,■) - q^fai-i), v qa .{ Vj ))- 

Here we denote the acute angle between the link q aj+1 (Vj+i) — Qvj (Vj) an d the 
inward unit normal v qa . +i ^ j+l) by Z( 9(T . +1 (^ + i)- ?(7 .(^), The condi- 

tion that -jf-L = is thus that the 2-link defined by the triplet (q^^ (<Pj-i, Qaj (Vj)) 1o- j+1 (Vj+i) 
is Sne// at </?j, i.e. satisfies the law of equal angles at this point. A smooth 
critical point of L on T M is thus the same as an M-link Snell polygon. Note 
that Snell polygons include polygons which have links inside the obstacles. 
They are known as 'ghost orbits' in the physics literature, since they cross the 
boundary between the interior and exterior domains. One may also include 
singular critical points of L a which correspond to polygons in which an edge 
is collapsed. 



2.2 Length spectrum 

We will denote by Lsp(Q c ), resp. Lsp(Q), the lengths of periodic billiard 
trajectories in f2 c , resp. Q. We do not include ghost orbits. By periodic 
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trajectory, we mean closed bicharacteristic curve; for the definition see [GM, 
PS]. Special periodic trajectories are m-link periodic reflecting ray, i.e. m- 
link Snell polygons. Their lengths are isolated in the length spectrum. 

When fl c is the complement of a pair of disjoint, strictly convex obstacles, 
the only exterior transversal reflecting rays are iterates Y of the bouncing 
ball orbit 7 between the components, i.e. the line segment between the points 
a± G 0±. When O is not assumed convex, there may be additional periodic 
reflecting rays and the shortest periodic orbit of Q c need not be the invariant 
bouncing ball orbit. Recall that we assume throughout that L 7 is isolated in 
Lsp(tt) ULsp(tt c ). 

It is part of the folklore of the Balian-Bloch approach that contributions 
to the trace formula due to ghost orbits are exactly cancelled by other ghost 
orbits. We will not prove that this cancellation occurs since the rigorous 
Poisson relation already implies that only Y contribute to the wave trace. 
Additionally, there are so-called diffraction orbits, which have points of tan- 
gential intersection with the boundary and may glide along the boundary. 
They are not critical points of the length functionals, and correspond not to 
smooth singularities in the wave trace, but to analytic singularities. 



3 Poisson relation for the exterior problem 

We now recall the Poisson relation for exterior domains in dimension 2. 

The exterior Green's kernel G^(k + ir, x, y) G V(VL C x Vl c ) with boundary 
condition B = D, or iV (Dirichlet or Neumann) is the kernel of the exterior 
resolvent R%(k+ir) = -(A%° + (k+ir) 2 )~ 1 : H s (tt c ) -> H s+2 (tt c ). Here, k G 
R% and r > 0. Its kernel (the exterior Green's function) may be characterized 
as the unique solution of the boundary problem: 



-(Ag c + (k + iT) 2 )G^{k + ir, x, y) = 5 y (x), (x, y G O) 

BG%(k + iT,x,y) =0, xedtt 

dG%(k+iT,x,y) \r*r (u ,. n A\ 

il k+tT)G R (k+tT, x, y) = o[ -), as r — > 00. 

or r 



(6) 



Here, the boundary operator could be either Bu = u\qq c (D) or Bu = d u u\g^c 
(N). We use a similar notation for the interior resolvent and boundary con- 
ditions. The boundary conditions D and iV are in a certain sense comple- 



11 



mentary, and we will write B' for the complementary boundary condition to 
B. 

We are interested in the regularized distribution trace of the combined op- 
erator Rg C (k+iT)Q)R B ,(k+iT) , more usually regarded as a distribution trace 
of the combined wave group E^f(t) © E^,(t), where E^it) = cost a/ — Ag c ). 
Let p E C$°((R) + ) and define 



R%{k + ir) ®R%,{k + ir) 

(7) 

p(k - p) (p + ir)R% c (p + it) © R%, (p + %T)dp. 



I 

Jr% 



From the resolvent identity (e.g.) 

1 f°° 

+ it) = — !— / e^ +iT ^E2\t)dt, 

fJ> + IT JO 



it follows that 



roc 

R%(k + it) © R%,(k + it) = / p{t)e* k+iT *E% c {t)®E%{t)dt. (8) 

Jo 

3.1 Poisson relation and scattering phase 

The following is the Birman-Krein formula for the relative trace of the wave 
group of an exterior problem, as generalized by Lax-Phillips and Bardos- 
Guillot-Ralston (see [BGR], Theoreme 3) : Let tp £ C^(R%). Then: 

Tr [ <p(t)[E%(t) © E%{t) - E (t)]dt = ±- [ <p(\)[ds B (\) + dN B ,(\)]. 
Jr% ^ Jr% 

Putting (p(t) = p(t)e^ fe+rr )* and rewriting in terms of the resolvent, we get: 



Tr[R%(k + it) © R% p (k + it) - R 0p (k + it)} 

= 7T [ p{k + iT-\)[ds B {\) + dN B ,{\)}. (9) 
^ Jr% 

Here, s B (X) = logdet Sb(A) where S B (E) is the scattering operator and N B > 
is the interior Weyl counting function with the indicated boundary condtions. 
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The singular support of the regularized wave trace Tr[E^f(t) — E (t)] 
is contained in the set of lengths of exterior periodic generalized billiard 
trajectories. That is, by [BGR], Theoreme 5, we have 

singsuppTr[E%(t) - E (t)} C Lsp{Q c ). 

Similarly for the interior wave group with Lsp(Q) replacing Lsp(Q c ). When 
L 7 is the length of a non-degenerate periodic reflecting ray 7, and when L 7 is 
not the length of any other generalized periodic orbit, then Tr[E^f (t) — E (t)} 
is a Lagrangian distribution in the interval (L 7 — e, L 7 + e) for sufficiently 
small e, hence Tr[E^(k + ir) — E p0 (t)} has a complete asymptotic expansion 
in powers of fc -1 . Let us recall the statement in the exterior case (see [GM], 
Theorem 1, and [PS] Theorem 6.3.1 for the interior case, and [BGR], §6 for 
the exterior case). Let 7 be a non-degenerate billiard trajectory whose length 
L-y is isolated and of multiplicity one in Lsp(Q). Let be a sufficiently small 
conic neighborhood of (M) + 7 and let x be a microlocal cutoff to T L . Then 
for t near L 7 , the trace of the wave group has the singularity expansion 

Tr X [E% c (t)-E (t)} ~ a 7 (t-L 7 +iO)- 1 +a 70 log(t-L 7 +iO) 

00 (10) 
+ J^a 7fc (t-L 7 +i0) fc log(t-L 7 +i0) 

fc=i 

where the coefficients a 7 fc are calculated by the stationary phase method from 
a Lagrangian parametrix. In the interior case, one of course omits the term 
E (t). 

The result may be re-stated as follows: Let p G C^°(L 7 — e, L 1 + e), equal 
to one on (L 7 — e/2, L 7 + e/2) and with no other lengths in its support. Then 
the interior trace TrR^ p (k + it) and the exterior trace Tr[R^ p (k + it) — 
Ro P {k + ir)] admit complete asymptotic expansions of the form 



Tr[R%(k + it) - Ro p (k + it)} ~ ^ 5 7;i A;- 



3=0 



TrR%(k + ir) ~ e l(fc+ir)rL - ^ 5 7;J A;- 



(11) 



whose coefficients B T j are canonically related to the wave invariants a 7; j 
of periodic (internal, resp. external) billiard orbits. We have removed the 
cutoff operator x since there are no singularities of the trace at t = L 7 in the 
microsupport of / — x- The coefficients depend on the choice of boundary 
condition but we do not indicate this in the notation. 
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3.2 Poisson relation and resonances 



We now recall the Poisson relation in dimension 2. Let {Xj} denote the 
resonances of A§ c , i.e. the poles of the analytic continuation of R^{k + it) 
from {r > 0} to the logarithmic plane A. Let 9 denote a (small) angle, and 
let Ag denote the conic neighborhood of the real axis in A consisting of points 
with angular coordinates in (— 6, 9). Further, let (p G Cq°((1R) + ) and let E (t) 
denote the free even wave kernel. The global Poisson formula (asserts the 
following: 

Proposition 3.1. (see Zworski [Zw], Theorem 1) 

Tr( [ m(E2 C (t) - E (t))dt = J2 m ( X M\) + "»(0M0) 

I nil * ' 

d A 3 eA e 

/°° dsn f°° 

i>{X) v {X)^-dx + J <p(t)is e ^(t)dt, 

with iy e ^eC°°(R%\{0}), dfu e ^ = 0(t~ N ), |i| ->oo, 

where sb denotes the scattering phase of A§ c , and where m(A) denotes the 
multiplicity of the resonance X. Also, tp G Cq°(R^ ) ) is a cutoff which equals 
1 for t near 0. 

(We note that the sum over eigenvalues term in the formula stated in 
[Zw] ( Theorem 1) is absent in the case of an exterior domain.) 

Substituting <p(t) = p(t)e^ k+iT ^\ with p G C£°(( R ) + ) with suppp suffi- 
ciently close to rL for some r, into (3.1) and using (8), we obtain: 

Tr[R%(k + ir) -R 0p (k + ir)} 

= miX^piX + k + ir) +0(|fc|- M ), (k^oo). ( 12 ) 

In the remainder estimate we use that 

m(0)p(k+tr), [ '^(X)p(X+k+ir)^dX, (u 6 ^, p{t)e^ kM ^) = 0(Ar°°). 
Jo dX 

In the case of the first two terms, this follows immediately from the fact that 
(fi(X) = p{k + it — X) and that both terms are integrals over a compact set 
of values of A. For the third term we additionally use that p vanishes in a 
neighborhood of t — 0; integration by parts and the estimate on on uq^ then 
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gives the rapid decay in k of the term. It follows that the asymptotics of the 
regularized resolvent trace Tr[R%(k + it) — Ro p (k + it)] are an invariant of 
the resonance poles. 

The following well-known proposition, explained to the author by M. 
Zworski, implies more: only poles in a logarithmic neighborhood {I^Ajl > 
iVlog \ Xj\} of the real axis cause singularities in the wave trace: 

Proposition 3.2. For any N > 0, there exists k(N) with k(N) — »■ oo as 
N — > oo ; such that 

J2 e C k(N \R% - {0}). 

= | QfA^- 1 <p| A^- 1 , 1 Of A^- 1 >AT log | | 

Proof. We have 



e iXjt \ 



A j : 1 9 Aj | < p | A j | , | A j | > JV log | A j | 

-Nt 



Aj : | S3f Aj | <p| Aj | , | S3f Aj- 1 > JV| Aj | 

< / r k ~ tN dN{r) <C I r^-^dr. 



(13) 



In the last line we used a polynomial bound on the number of resonance 
poles (cf. [Zw] (2.1) for the estimate). Thus, the portion of the wave trace 
coming from poles outside of sufficiently large logarithmic neighborhoods of 
the real axis is as smooth as desired. □ 



4 Neumann expansion and Poisson relation 

In the Balian-Bloch approach to the Poisson relation, one uses the classical 
theory of layer potentials ([Til] §7. 11) to express G^(k + ir,x,y) in terms 
of the 'free' Green's function G (k + iT, x, y) of (M) 2 , i.e. the kernel of the free 
resolvent (A — k + ir 2 )' 1 of the Laplacian A on (R) 2 . A classical reference 
for this approach is the paper of Pleijel [P]. This approach is also used in 
[M] to make a reduction to the boundary in studying resonance poles and in 
the duality result of (3). We review the necessary background and notation. 
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4.1 Potential theory 



We recall that the free Green's function in dimension two is given by: G (k + 
ir, x, y) = H^\k + ir\x — y\), where H^\z) is the Hankel function of index 
0. We recall its small and large distance asymptotics ([TI], Chapter 3, §6): 



"2^ ln(r) 
iF(v) ,2 



TT 



as r — > 0, if v = 



Y as r -> 0, if v > 



(13) 



3 j(fe+«Tr— viv/2— it/ 4) 
r l/2 



1 + O (r-fg) j as |k + ir|r -> oo. 



Here it is assumed that r > 0. Thus, Go(k + ir, x, y) has two kinds of asymp- 
totics: a semiclassical asymptotics for large (k + ir)\x — y\ and a singularity 
asymptotics along the diagonal. It is the first kind of asymptotics which 
accounts for the connection to billiard trajectories. The singularity along the 
diagonal will have to be regularized. 



4.1.1 Layer potentials and boundary integral operators 

The double layer potential is the operator 

V£(k + ir)f(x) = [ ^-G Q (k + it, x, q)f(q)ds(q), (14) 
Jan ov y 

from H s (d£l) — > H^ 1 ^ 2 ^), where ds(q) is the arc-length measure on dQ, 
where v is the interior unit normal to Q, and where d v — v • V. It induces 
the boundary operator 

N(k + ir)f(q) = 2 / j-G (k + ir, q, q')f(q')ds(q') (15) 
Jon ov y 

which map H s (dtt) -> H s+1 (dQ). 
By the explicit formula we have: 

^N(k + ir,q(ip),q(ip')) = d Uy G (^,q(ip),q(ip')) 

= _ (k + iT)H?\k + iT\q{ip)-q{v')\) 
x cosZ(g(y?) - q(tp'),v q{ip) ). 
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We will need the 'jump' formula of potential theory ( [TI] Chapter 7. 11, 
(H-7)), 

(V£(k + ir)u) ± (q) = ^(N(k + ir) =F I)u(q), (16) 

where 

f + (x) := \im x ^ q:Xe0 f(x), 

(17) 

f-(x) = \im x ^ qjXe nf(x). 

Note that the sign of each term in (16) depends on the choice of the 
interior/exterior: in N(k + ir) it depends on the choice of interior /exterior 
unit normal and in the identity operator term it depends on whether the 
limit is taken from the interior or exterior. In the interior case we have 
(V£(k + ir)u)±(q) = l(N(k + ir) ± I)u where N{k + ir) is defined using the 
interior normal and where ± have the same meanings as in the interior case. 

4.2 Interior/Exterior duality 

In [Zl, Z2] we used the classical reduction of the Dirichlet problem to the 
boundary to study the Dirichlet resolvent (as in [BB1, BB2]). We now use a 
similar method to obtain a convenient formula for the Fredholm determinant, 
det((J + N{k))). Similar ideas can be found in [THS, THS2, AG, GP]. 

The following formula is sometimes referred to as interior/ 
exterior duality. Combined with the Birman-Krein formula (9) it expresses 
the growth rate of the inside plus outside spectra in terms of the determinant 
of a boundary integral operator: 

Proposition 4.1. For any r > 0, the operator (I + N(k + ir)) is of trace 
class and has a well-defined Frehdolm determinant, and we have: 

Tr (Rf [R%(k + ir) © R% - R 0p (k + ir)] 

f d 
= / p(k - A) — logdet(J + iV(A + ir)d\, . 
J R n dX 

where det(I + N(\ + ir) is the Fredholm determinant. 

Proof. We first argue formally. The interior /exterior resolvent kernels can be 
constructed in the classical way, as follows (see e.g. [P]). For simplicity, we 
consider the interior Dirichlet resolvent, but the construction for the exterior 
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resolvent or with a change to Neumann boundary conditions is almost the 
same. We have: 

R% c (k+ir) = lnciRoik+i^-Viik+i^il+Nik+iT^-^Roik+iT^lnc, 

(R%(k+ir)) tr = InlRoik+i^-Viik+iT^I+Nik+iT^-^R^k+iT)}!^ 

(18) 

where 7 denotes the restriction to the boundary taken from within the rel- 
evant domain. When we take the regularized trace, we first subtract from 
R^(k + ir) © R% c (k + ir) the free operator Ro(k + ir), which removes the 
first terms on the right. We then cycle the factor ^R^k + ir) from the right 
to the left side, obtaining an operator on dVL. The inside and outside terms 
add up to the kernel 

/ G (k + iT,q',y)d I/q G (k + iT,y,q)dy. (19) 

This indeed is why the interior Dirichlet and exterior Neumann problems 
were combined and explains the sense in which they are complementary. We 
now claim that the resulting kernel equals ^^N{k + ir, q', q). To prove this 
we note that 

-^rRo(X) = 2A(A + A 2 )- 2 = 2Ait> (A) 2 . 
of A 

Hence, 

^iV(A) = A 7 ^ o(A ) 7 = 2\ 1 d Vq R Q {\) o J R (A)7 

which is precisely the kernel (19). Hence the right side of (18) equals log(/+ 
N(k + ir)). The final formula follows by putting k = A, integrating against 
p(k — A) (A + ir) and noting the cancellation of the second factor. 

To justify the formal manipulations, we need to show that the relevant 
traces and determinant are well defined. This will be done in a sequence of 
Lemmas. 

Lemma 4.2. For any r, N{k + ir) G Xi(dQ), the ideal of trace class oper- 
ators on L 2 (dQ). Hence, det(J + N(k + ir)) is well-defined as a Fredholm 
determinant. 

Proof. In [Zl] it is proved that N(k + ir) e ^~ 2 (d£l) and this immmediately 
implies that it is of trace class. In fact, N(k+ir, q, q') has just a \q—q'\ log \q — 
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q'\ singularity on the diagonal. That it is trace class also follows from the 
Hille-Tamarkin theorem (see [Zl] for references). 

It is then a classical remark that the Fredholm determinant det(J + N(k + 
ir)) is well-defined, in fact | det(7 + N(k + < exp(| |iV(fc + «t) | |i), where 
|| ■ ||i is the trace norm. See [S], Lemma 3.3. 



Remark 4.3. 1 Note that the usual statement (cf. [Til]) is that N(k + it) 

is of order —1. However, the principal symbol vanishes in dimension 2. 

2 We recall (see [S], Theorem 3.10) that the Fredholm determinant is given 



det(7 + N(k + ir)) = — / det[N(k + it, q h qj)]i<i,j< n ds{q 1 ) ■ ■ ■ ds(q n ). 



Again combining classically known facts, we have: 

Lemma 4.4. For any r > 0, logdet N(k+ir) is well-defined and differentiable 
in k, (I + N(k + ir^^N'ik + ir) is of trace class and we have: 

log det N(k + it) = Tr 9n (l + N(k + iT))~ l N\k + it). 

Proof. First we note that (I + N(k + it)) is invertible on L 2 (dVl) if r > 0, so 
det (J + N(k + it)) ^ 0. Hence its logarithm is well-defined. Differentiability 
of logdet(J + N(k + it)) and the formula for the derivative follows from the 
general fact that det(J + A) is Frechet differentiable on Ti with derivative 
equal to (/ + A)' 1 if —1 cr(A) (here, c(A) is its spectrum; see [S], Corollary 



The singularity on the diagonal of N(k + iT, q, q') is independent of k + tT, 
so the derivative is at least as regular. Hence, N'(k + iT) e X\. The statement 
that (J + N(k + iT))~ l N\k + it) e Z x follows from the fact that (I + N(k + 
it)' 1 E ^°(dtt). 



□ 



by: 




5.2). 



□ 



Lemma 4.5. The operators 



(?) l n J R n p(k - A)W(A+*T)(/+A^(A+*r))-VRo(A+*r)]dAi n , 
(ii) l Q c j Rn p(k - A)W(A+2r)(/+iV(A+2r))- 1 7i? (A+2r)]rfAl nc , 
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are of trace class and the sum of their traces equals Tra(I+N(k+ir)) 1 N'(k+ 
ir). 



Proof. We expand (J+A^A+ir))- 1 = I - N(X)(I+N(X+ir))- 1 . The / term 
gives us T>£(X-Ht)^R (X-Ht)]. We break up each of V£(X+ir) and ^^(X+ir)] 
into regular (continuous) and singular parts using the small z expansion 
(13) of the Hankel functions Hq(z), Hl(z). The singular parts of the kernels 
are independent of (k + ir). Hence the product of the singular parts gets 
multiplied by p(0) = upon integration against p(k — A). Therefore only 
the regular parts survive the convolution, and these are of trace class. The 
remaining term is a bounded operator composed with trace class operator 
for each A, so it is of trace class for each A. 

We may then cycle the factor of V£(X + it) (or each of its regular and 
singular parts) to the right side of the traces. After doing so, we reassemble 
the regular and singular parts into r yR ( y X + ir)]V£(X+iT)and we reassemble 
I - N(X)(I + N(X + ir)y 1 into (/ + JV(A + zt)) -1 . The calculation at the 
beginning of the proof then gives the stated formula. 



Corollary 4.6. Suppose that L 7 is the only length in the support of p. 
Then, 



where B r j are the wave invariants of j in (11). 

To prove Theorems (1.1) and (1.4), it thus suffices to determine Vt from the 
integrals in Corollary (4.6). Since we have combined the interior/exterior, we 
emphasize that our inverse result only assumes knowledge of the resonances: 
Indeed, suppose the resonance poles of Q c are known. Then by Proposition 
(3.1), the asymptotics of the exterior resolvent trace in terms of exterior 
periodic orbits are known. But by Proposition (4.1), these are the same as 
the asymptotics of the integrals in Corollary (4.6). 



□ 



This concludes the proof of Proposition 4.1. 



□ 
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5 Trace asymptotics 

We now explain how to use Corollary (4.6) to calculate the coefficients Byj. 
We write 

logdet(J + N(X + it)) = Tr dn (I + N(X + ir))- l N'(k + tr), 

dX 

and then expand (J + N{k-\-iT))~ 1 in a finite geometric series plus remainder: 
(I+NiX+ir))- 1 = ^(-1) M iV(A) M + (-l) Mo+1 iV(A) Mo+1 (/+iV(A+*r))- 1 . 

M=0 

(20) 

We now argue that for each order k J in the trace expansion of Corollary 
(4.6) there exists M (J) such that 

(0 £m°=o(-1) M 7> k% P( k - A ) N(X) M N'(k + ir )dX 

(ii) Tr j Rfl p{k - X)N{X) Mo+1 {I+N{X+iT))' 1 N'{k + ir)dX = 0(k~ J ~ l ). 

(21) 

We will sketch the proof of (i) in some detail since it explains how to 
calculate the coefficients B 7 rj. For a discussion of the remainder (in an 
analogous but not identical calculation) we refer to [Z2]. 

To simplify the notation, we integrate by parts in (i) to throw the deriva- 
tive onto p, and then the issue is to analyse the traces: 

Tr [ p(k- X)N(X + ir) M dX 
Jr% 

= I I p(k- X)[Uf =l N(X+ir,q vqj+1 )ds(q 3 )]dX (where q M+1 = i$) 
Jr% J(dn) M 

Since 

N(k + ir) : L 2 {0+) © L 2 (0_) -> L 2 (0 + ) ® L 2 {0-) 
we write it as 



N(k) 



N ++ (k) N + _(k) 
N_ + (k) N^_(k) 
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Then (23) is a sum of terms m}->{±i} ^m,p wrtn 

Im,p= [ I p(A;-A)[n^ 1 iV (Tji(7j+1 (A+.r,g (Tj (^),g^ 1 (^ + i)^]^. (24) 

Here and hereafter, indices are understand modulo M. Using the asymp- 
totics of the free Green's function in (13)-(27), the crth term is formally an 
oscillatory integral with phases given by the length functional (4). We would 
like to calculate the trace asymptotically by the stationary phase method. 

As discussed extensively in [Zl], we cannot immediately do so because 
of the singularities along the diagonals of the integrand. We therefore first 
need to de-singularize the integrals. One of the main advantages of the 
present inside/outside reduction to the boundary is that it simplifies the 
regularization procedure by eliminating the integral over Q in [Zl]. We briefly 
outline the method and refer the reader to [Zl] for further details. 



5.1 Boundary integral operators as quantized billiard 
maps 

The operators N_ + (k + ir),N + _(k + it) are semiclassical Fourier integral 
operators with phases \q+(<p) — q-((p')\. The phase (and amplitude) is non- 
singular since the boundary components are disjoint. The phase in the H — 
case parametrizes the graph of the following multi-valued billiard map from 
B*0 + to B*0-\ given (q+,v+) G T q O + with \v\ < 1, add a multiple of 
the unit outward normal to turn v into an outward unit vector and proceed 
along the straight line in that direction and let g_ be a point of intersection 
of this line with O-. Let V- be the tangential projection of the terminal 
velocity vector at g_ , and put (3 + -(q + , v+) = i>_). We have (deliberately) 
described in an ambiguous way: when q_ is the first intersection point, 
we have the usual exterior billiard map; but the phase actually parametrizes 
the canonical relation which includes intersection points which occur after 
the line enters the interior of O-. As mentioned above, these are 'ghost 
orbits' of the billiard flow and they cancel out of the trace formula and do 
not contribute to the wave invariants we are calculating. 

The operators N ++ (k + it) and (k + it) are more complicated, and 

also, as it turns out, less important. Each is a combination of a homogeneous 
pseudodifferential operator of order —1 with singularity on the diagonal and 
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a Fourier integral operator of order which 'quantizes' the billiard map of 
the interiors of Q and Q c . They are exactly the kind of operators discussed 
extensively in [Zl, HZ]. There is a further transition region between these 
two regimes in which it behaves like an Airy operator, but this region will 
not be important for our problem. 

To separate the two basic regions (tangential and transverse) of iV ++ and 
N , we introduce a cutoff of the form x(/c 1_<5 (y2 — if')) in terms of the arc- 
length coordinate tp used in the parameterization q(ip). Here, \ a smooth 
bump function which is supported on [—1,1] and equals 1 on [—1/2,1/2]. 
We then write (dropping the subscripts) 

N(k + it) = N (k + it) + N x (k + ir), 

< N (k + it, q(<p), q(<p')) = xik 1 - 5 ^ - <p'))N(k + it, q(<p), q(<p')), 

k N 1 (k + it, q( V ), q{y>)) = (1 - xik 1 - 5 ^ - tp'))N(k + it, q{ip), q(<pf)). 

(25) 

Each term is a (non-standard) Fourier integral operator with amplitudes 
belonging to the following symbol class: We denote by Sg(T m ) the class of 
symbols a(k, <pi, . . . , ip m ) satisfying: 

\(h- 1 D v ) a a(k,<p)\<C a \k\ p - 5 M, (\k\>l). (26) 

This follows from the asymptotics of Hankel functions (13), which may be 
described as follows: there exist amplitudes a , a± satisfying 

(1 — x(k l ~ s z))ai({k + %t)z) G S$(R%), 

(l-x(k 1 - 5 z))a ((k + l T)z)eS; 1/2 (R n D ) 

and such that 

Ui) H^({k + %t)z) = e^ h+lT >a ((k + it)z), 

I (27) 

((H) (k + iT)H^\(k + iT) Z ) = (k + iT)h- l({mT)z ai ((k + iT)z), 

The iV term has a singularity on the diagonal of a pseudodifferntial 
operator of order —2. The Ni term is manifestly an oscillatory integral 
operator of order with phase \q(ip) — q((f')\- As discussed in the — h case, 
it generates the interior billiard map (interpreted in the multivalued sense of 
the — h case). 
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5.2 Microlocalization to 7 



So far, our considerations have been global. But our goal is to calculate wave 
invariants at a bouncing ball orbit, and for this purpose we can microlocal- 
ize the trace to a neighborhood of 7. This obvious sounding statement is 
somewhat problematic in the present approach since the relevant operators 
are not standard pseudodifferential or Fourier integral operators. 

The billiard map is a cross section of the billiard flow, and in this 
cross section 7 corresponds to the periodic points (0, 0) G B*0 + and (0, 0) e 
B*(0-) of period 2 for Here, we choose the parametrizations of 0± 

such that if = at each endpoint of the bouncing ball orbit; it is normal 
to 0±, so its tangential part is the vector (0,0). To microlocalize to this 
orbit we introduce a (block diagonal) semiclassical pseudodifferential cutoff 
operator 



with complete symbol xifiV) supported in V e := {(</?, 77) : \ip\, \r)\ < e}. We 
also write 

p*4r logdet(J + N(k + ir)) = Trp * (I + N(k + it))' 1 o 4r N ( k + ir )- 
dX dk 

Recall that to test whether (L, r) e WF(u) (wave front set) for u(t) £ 
S'(R%), one checks whether F(pu)(k) ~ for all p E C^°(R^)) satisfying 
p = 1 in some interval (L — e, L + e). Here, Tu = u denotes the Fourier 
transform. Equivalently, one checks whether p*u ~ 0. We now show that in 
calculating wave invariants we can microlocalize our non-standard operators 
to 7 in the standard way. 

As a preliminary step, we prove the simpler relation 



p * Tr[V(k + ir) (/ + N(k))- ltr S(k + ir) tr o (1 - x 7 (Jfe)] ~ 0. (28) 



where x~,(k) is a semiclassical cutoff in Q c to a microlocal neighborhood of 7 
which has the form 



near the boundary. Here, (r, y) are Fermi normal coordinates in Q c near the 
endponts of 7 , with r the distance to the boundary. Thus, differentiations 
are only in tangential directions. Note that (28) is equivalent to saying that 




Xf(ty = X(r,y,k l D y ) 



Tr{R° c D -l nc R 0p ]o(I-x y (k))~Q, 



(29) 
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In view of (18), (29)-(28) follow from standard wave front set (WF) 
considerations for the exterior even wave group E^ c (t,x,y), and from our 
assumption that Y is the only exterior periodic orbit of its length. Recall 
that WF(E% C (t, x, y)) is the space-time graph 

T = {(t,T,x^,x>,0 £T*(R%xn c xW) : t = G\x, = (x\ £')}, 

(30) 

of the generalized billiard flow in the exterior. The only unusual feature is 
our particular choice of cutoff operator to microlocalize to 7. To verify that 
it behaves correctly like a microlocal cutoff to 7 we write 

POD 

R%(k)x 1 (k) = / p{t)E^\t)^y,\D t \- l D y )e i{ - k+i ^dt. (31) 
Jo 

Now in Fermi coordinates (r, pdr, y, rjdy), (R) + 7 is ray in the direction of dr. 
Since the space-time graph of the cotangent bundle along 7 may be described 
in normal coordinates as a neighborhood of 

{(t,T,t,T,0,0)} 

a conic neighborhood may be described in these coordinates by \y \ < e, \t}/t\ < 
e. This is precisely the set to which x 7 (r, y, \D t \~ l D y ) microlocalizes. Thus, 
(29) and hence (28) are correct. 

We now show that this empty WF relation remains true after moving 
the cutoff to the left of S£(k) tr as a cutoff pseudodifferential operator on 
the boundary to a neighborhood of the billiard map orbit. As mentioned 
above, this would be straightforward if T>£(k) and the other operators in the 
composition were standard FIO's. However in addition to their semiclassical 
canonical relations they have singularities on the diagonal and both con- 
tribute to their WF. We will shorten the notation for the cutoff operator to 

Xo- 

Lemma 5.1. We have: 

Trp * (/ + iV(A; + ir))' 1 o ^N(k + it) 

~ Trp * (I + N(k + ir))' 1 o Xo (k) o ±N(k + it). 
Proof. We first express the various layer potentials in terms of the wave 
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group: 

' S£{k + ir) tr = 7 Go(A; + it) = J °° e^ k+iT ^E (t)dt, 
< l Q c(V£(k + it) o (/ + N(k + ir))- 1 ^, q) = P%°(k + it, x, q) 

k = lv2G^(k + ir, x, q) = f °° e^ k+i ^ ll/2 E^(t)dt, 

where Pjf is the Neumann Poisson kernel of the exterior domain and where 
j u u = dviilonc^u = u\gnc. We have subscripted the restriction operators 
to clarify which variables they operate on. The composition T>£{k) o (7 + 
N(k))^ 1 o (1 — Xl (k)) ° S£(k) tr may be written (with the relevant value of k) 

as 

rOC nt 

/ e it(k + ir) { lv2 E%(t-s)o(I- Xl )(y,\D t \- 1 D y ) 1 oE (s)ds}dt. (32) 
Jo Jo 

We therefore have: 
p * S£{k) o (/ + iV(A;))- ltr o (1 - Xl (k)) o V£(k) tr 

= Jo°°P(t)[Jo^E% c (t -s)o(I- Xl )(y, IAI-^^T.i o E (s)ds]e i ^ t dt. 

(33) 

Our goal is to show that 

WF[Tr [ lu2 E^(t-s)o(I- Xl )(y,\D t \- 1 D y ) 1 oE (s)ds}n[rL 1 -e,rL 1 +e} = 0. 
Jo 

(34) 

This follows as long as 

V(s,t) = Tr lv2 E%(t - s) o (/ - Xj )(y, lAl^A/h o E (s) 

is a smooth function for t G (rL 7 — e, rL 7 + e) and for s G (0, rL 7 + e). To 
prove that it is smooth, we observe that the singular support consists of (s, t) 
such that there exists a closed billiard orbit of length t outside a phase space 
neighborhood of Y which consists of a straight line segment of length s from 
a point x G Q c to a boundary point q, followed by a generalized billiard orbit 
of length t — s from q back to x. By our assumption on the length spectrum, 
the only possible orbit with length t G (rL 7 — e, rL 1 + e) is 7 r of length r; 
but the cutoff has removed this orbit. This completes the proof. 

□ 
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5.3 Regularization 

In the preceding section we explained that the trace can be microlocalized 
to 7. However, it is not true that the operator N can be microlcalized to 7, 
that is, we cannot replace N by x-y^X-y- This is due to the singularities on 
the diagonal. We now explain how to regularize them. 
When we expand 

(N ++0 + N ++l + N + _ + N__ + N__, + iV_+) M , (35) 

we get a sum of 'strings' of factors. The factors of iV ++ o or iV forbid the 

immediate application of the stationary phase method due to the singularities 
of its phase and amplitude on the diagonal. We now explain how to remove 
such factors and regularize the term except for presence of N+ +0 or N M _ , 
which has to be handled separately. We will only explain how to remove one 
factor of iV ++0 next to a smooth factor such as N ++1 . In general, one just 
repeats the procedure one pair at a time. 

5.4 The compositions Nq o Ni 

We begin by characterizing the composition N o Ni. Essentially the same 
result holds for Ni o N by the same argument. 
The composed kernel equals 

N oN 1 (k + it, <p u y? 2 ) := (k + irf J T - <p 3 )){l - x{^ 1+& ^ 2 - <p 3 ))) 

H ( i \{kfi + ir)\q(ip 3 ) - q((fi)\) cosZ(g(v? 3 ) - qM, ^(^)) 
H[ 1] ((kfi + ir)\q((p 3 ) - qM\) cos Z.{q{ip 2 ) - qM, ^ 2 ))^3- 

(36) 

The next proposition shows that these compositions are semiclassical 
Fourier integral kernels. Thus, the singularity on the diagonal is 'regular- 
ized'. The crucial point is that the composition with the singular factor iVo 
decreased the order by one unit. We state in a somewhat simplified form and 
refer to [Zl] for a more detailed statement. 

Proposition 5.2. N oNioxo(k+iT, <pi, (£2)) defines a semiclassical Fourier 
integral operator on dfl of order —1 associated to the billiard map. 
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Proof. We only sketch the proof and refer to [Zl] for further details. 

The proof is an explicit calculation. Following [AG], we change variables 

and then to change variables 1$ ' — > w, with: 

{ItfM ¥>i > V?3 f — 1?) ^>o 

-k(^s) - g(^i)|, ^1 < ^3 [ -#) -q(<pi)\, ^<o 

(37) 

In other words, we change from the intrinsic distance along dQ to chordal 
distance. The purpose of this change of variables is to simplify the difficult 
factor in (36): 

+ tr)\q( V3 ) - q(<pi)\) -> + 

We then substitute the asymptotic WKB formula for the 'easy' factor 
H^\{kn + ir)\q(ipi) — ^(^3 1 ) which is valid on the support of the cutoff. The 
key points are then that 

(ii) cos(Zg((/9 2 ) - <?(y?2 + v q{^ 2 )) ~^ \ u \K((p 2 , u), with K smooth in u; 

< {ii)H [ l\{k^ + ir)\q( Vl ) - q( V2 - -> e ik \^- q ^e ikua A{k + it, (p u <p 2 , u), 

where Ak is a symbol in k of order —1/2 and smooth in u. 
Here, 

a = sin $i i2 , with i?i j2 = Z(g(v?2) - ?(v?i)> ^2))- ( 38 ) 

It follows that the composed kernel (36) can be expressed in the form 
Ae i(k+ir)\q(<pi)-<P2)\ (f urt her composed with Xo) with 

A(k + ir, <p u y? 2 ) = JZ X{k l - 5 u){\ - xik 1 - 6 ^ -<pi- «)) 

(39) 

x \u\e ikau H[ 1] ((k + ir)\u\)G((k + ir),u, </?i, V2 )du, 

where G((k + ir),u, tpi, ip 2 ) is a symbol in k of order —1/2 and smooth in u, 
and where x(/c 1_<5 -u) = x(A; 1 ^ 5 ((/?i — y?3)). 
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We now change variables again, v! = ku (and then drop the prime), to 

get 

A(k + ir,<p 1: <p 2 ) = fc- 2 /^ x(fc"*«)(l-x(fc 1 "*(¥'2-¥'i-A;- 1 «)))He <o « 
• H[ 1 \b\u\)G(k + tr,l Vl , V2 )du, 

(40) 

with b = 1 + i{r/k). Since \u\ < k s on the support of the cutoff, we have 
||| < k~ 1+s . The Taylor expansion of G{k + ir, u, ip±, ip 2 ) at u — produces 
an asymptotic series in A;. The main point is that the change of variables 
introduced a factor of k~ 2 , which cancels the original factor of k 2 in (36). We 
are left with an amplitude of order —1/2, which defines a Fourier integral 
operator of order —1 with phase \q(<fi) — ^2) | as long as the integral (40) has 
the same order as G. 

After Taylor expanding, we end up with terms of the form 

/oo 
x(k- s u)\u\u n e iau H^\b\u\)du 
-00 

times simple functions of (ipi, ip 2 ). The integral may be expressed in the form 

JZo x(^ 5 «)e M «^ 1) (6| M |)^| a=sin „ li2 , b=(1+Wfe) . (41) 

Using the cosine transform of the Hankel function, one can explicitly evaluate 
(41) as J^r(l — a 2 )~ 3 / 2 , modulo lower order terms. After composition with 
the cutoff xo this factor is bounded above, so one gets that A has order —1/2. 

□ 

We iterate this proposition to deal with the full strings in (35). The 
only cases where this is not possible are the ones which have no factors of 
N + _(k), N_ + (k). However, they are still composed with the cutoff operator 
X 7 and a somewhat similar calculation shows that this composition produces 
a semiclassical pseudodifferential operator, whose order is decreased from 
by one unit of k for each factor of N , 

5.5 Final remarks and remainder estimate 

After removing the Nq factors from the terms (35) of (20, we end up with a 
finite sum of ordinary oscillatory integrals plus a complicated remainder. In 
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the next section, we will describe the result of applying stationary phase to 
the finite sum. We now make a few comments on the remainder, referring to 
[Zl] for further discussion. 

The estimate of the remainder is complicated because the operator norm 
| \N(k+ir) 1 1 does not decrease with increasing r. To obtain a small remainder 
we redefine r — > r log k. This changes the wave trace expansions by k~ CrL ~< 
and hence for R sufficiently large a remainder estimate of 0(k~ R ) is sufficient 
to separate a finite part of the wave trace from the remainder. We then 
estimate the remainder 



by applying the Schwarz inequality for the Hilbert-Schmidt inner product. 
We futher use standard estimates on the Poisson kernel to remove the factor 
of (J+A^(A+?r))~ 1 . Unfortunately, this estimate will also remove the cutoff 
p and replaces N M by N M N* M . If one regularizes these products as above, 
we find that the critical points correspond to closed circuits of M-links which 
begin at some point x, end at some point x' and then return to x by travers- 
ing the links in reverse order. The cutoffs xo 5 X7 prohibit a proliferation of 
small links (grazing rays) and force the links in critical paths to point in the 
direction of 7 and hence to be of length roughly ML 7 . The imaginary part 
ir log k of the semiclassical parameter then contributes a damping factor of 
e -rML 7 log k f or eac h ii n ]^ The links correspond to the Ni factors. Thus, for 
each string, we have one k~ l for each Nq factor and one e - TML i lo & k f or each 
N\ factor. For sufficiently large r these combine to give a factor of k~ R for 
any prescribed R for each term of (35). 

6 Calculating coefficients 

We now apply the stationary phase method directly to each regularized in- 
tegrals. Additionally, we have to prove that the remainder is as small as 
claimed. The result is the following explict formula for the wave invariants: 



Tr / p(k-\)N(X) Ma+l (I+N(\+iT)Y l Xl N'(k + iT)d\ 
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Proposition 6.1. We have: 

B YJ + B 7 -rj_i = b r , j:2j f 2j (0) + b r ^ 2j -if j ~\ti) + B YJ , with 



2r 

- {2(/lll)J 2-2cosa/2 + B^WW^O)}} 

q=l 



where B^rj depends only on the 2j — 2 Taylor polynomial of f at x = 0. 

We review the calculation of the coefficients from [Z2] and adapt it to the 
exterior case. We also take advantage of the reduction to the boundary to 
simplify the arguments. 

6.1 Setting things up 

We recall that, in a small strip T e (j) around a + a_, the boundary <9f2 consists 
of two components which are symmetric graphs over the x-axis. We write 
the graphs in the form y = ±/(x) near a±. 

As in (23), integrals over (Q fl T e (7)) M consist of 2 M terms, corresponding 
to a choice of an element a of 

{±} M := {a : Z M -> {±}}. 

The length functional in Cartesian coordinates for a given assignment a of 
signs is given by 

2r-l 

c a (x 1} . . .,x 2r ) = yj ( x j+i- x j) 2 +(f«u+i)( x i+i)-fvU)( x j)) 2 - ( 36 ) 

If we write out the integrals of (24) in Cartesian coordinates, we obtain 

p poo pco f 

J(-e,e) M JO J-oo I 

X \{ XpH - X p , fa{pH)( X pH) ~ f*( P )( X p))\) ^ 
(XpU— X p)~ fa(j^l)( x pH){f(T(pH)( X pH)~ fa(p){ x p)) 
yj (Xpa - X p ) 2 + {fa(pn) { X pH ) + fa(p) (XpH ) (x p ) + L) 2 

x x (l - A)e ifc(1 - A) *rfArftrfxi . . . dx M - 
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We now regularize the integrals (37) as above and obtain classical oscil- 
latory integrals in x — y coordinates. We also eliminate the (t, A) variables 
by stationary phase. The Hessian in these variables is easily seen to be non- 
degenerate, and the Hessian operator equals — ^q\- Since the amplitude de- 
pends on t only in the factor p(t), which is constant in a neighborhood of the 
critical point, only the zeroth order term in the Hessian operator survives the 
method of stationary phase. We may therefore eliminate the dsdp integrals 
and replace the amplitude and phase by their evaluations at p, = l,t — rL 1 . 
This simplifies the above to 

k~ us f he ^£ ikC±( - Xl >---> X2m) p{£±{xi, x 2m ))a(k, x 1 ,x 2 ,..., x 2r )dxi ■ ■ ■ dx 2r , 

(38) 

An examination of the amplitudes and of the regularization procedure 
leads to the following conclusions: 

• (i) The regularized integral I M , p {{k + it)) is negligeable as k — > oo 
unless M = 2r (where rL 1 is the unique length in the support of p), 
and where only strings of iV_ + and iV + _ are left upon regularization. 
Otherwise there are no critical points. 

• (ii) The amplitudes of these oscillatory integrals have the form 

oo 

^ ] n A n {x)i 

71=0 

where A n depends only on the first n + 2 derivatives of f a . 

• There is a special term, namely the initially regular one in which 2r 
factors of N_ + and iV +1 alternate, corresponding to 2r bounces of the 
bouncing ball orbit. The phase is simply C± and the amplitude is 

a± (k,X!, . . . ,x 2r ) = 

n^ 1 ai((/c+ir)/iy / (a; p -a;, >fl ) 2 + (/ (7± (p)(a; p )-/ (7±( ^ 1) (a; ?H _ 1 ) 2 ) 

(x p — Xp+l)f'( p )(Xp) — {fa±(p)(Xp) — fa±(p+l)(Xp+l)) 

X 

^/(Xp - X p+1 ) 2 + (fa ± (p)(Xp) ~ fa ± (p+l)(x p+1 )) 2 

(39) 

We only need the principal terms in the symbols a±, but to simplify 
notation we omit this step. 
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This amplitude and phase of the special term have certain key attributes 
which will be used extensively below 
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(i) In its dependence on /, the amplitude has the form 
A(x,y,fJ'). 

(ii) In its dependence on the variables x p ,the amplitude 
has the form: 

a°(k, xi,..., x 2r ) = U^ =1 Ap(x p , x p+1 ) (2p + 1 = 1) 
< (in) Dfr l) C^{{x p -x p+1 f + {f ev {x p ) m 

~ fe p+1 (Xp + l)) 2 )- 1/2 (fe P (x P ) 
~ fe P+ Ax P+ l))f£ j 

k (iv) Va°(k,x 1 ,...,x 2r ) x =o = 0. 

Here, the notation = means equivalence modulo terms with fewer derivatives 
of /±. The last statement follows from the values 

/(0) = L, /'(0) = x\ x=0 = 0, Z(A - B, u B ) = tt/2, (41) 

where A, B denote the endpoints. 

We now use this information to determine where the data f±(0), 
yft?- 1 ) g rg ^ appears in the stationary phase expansion for I k (a So , C±). 



6.2 Stationary phase expansion 

We now apply the method of stationary phase to these integrals and keep 
careful track of the dependence on the number of derivatives of f±. To analyze 
the jungle of terms in the stationary phase expansion, we use the diagram- 
matic approach (see e.g. [AG]). We give a brief review. 

Consider a general oscillatory integral = f, R \n a(x)e' lkS ^dx where a G 
C^°((lR) n ) and where S has a unique critical point in suppa at 0. Let us write 
H for the Hessian of S at 0. The stationary phase expansion takes the form: 



a_\ n/2 iir S gn{H)/i 
Z7r \ f e ikS(0) 2?h£ 

k) yJ\detW 



where Z«=E^{ E J^} 



3=0 ^ (r,£): Xr / 
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Here, the sum runs over the set Gyj of labelled graphs (T, £) with V 
closed vertices of valency > 3 (each corresponding to the phase), with one 
open vertex (corresponding to the amplitude), and with / edges. Further, 
the graph T' is defined to be T minus the open vertex, and xv = V — I equals 
its Euler characteristic. We note that there are only finitely many graphs for 
each x because the valency condition forces I > 3/2V. Thus, V < 2j, / < 3j. 

The function £ 'labels' each end of each edge of T with an index j e 
{1, . . . , n}. Also, S(T) denotes the order of the automorphism group of T, and 
Ze(r) denotes the 'Feynman amplitude' associated to (r, £). By definition, 
h(T) is obtained by the following rule: To each edge with end labels j, k one 
assigns a factor of ^rhP k where if -1 = (h jk ). To each closed vertex one assigns 

d v S( 0) 

a factor of ik dxh „; dxiv where v is the valency of the vertex and i\ . . . , i v at 
the index lables of the edge ends incident on the vertex. To the open vertex, 
one assigns the factor dxil "^v ■, where v is its valence. Then Ie(T) is the 
product of all these factors. To the empty graph one assigns the amplitude 
1. In summing over (T,£) with a fixed graph T, one sums the product of all 
the factors as the indices run over {1, . . . , n}. 

6.2.1 The special term: The data /± 5 (0) 

We first claim that /i 2i) (0) appears first in the k j+1 term. This is because 
any labelled graph (r, £) for which h(T) contains the factor f± j \o) must 
have a closed vertex of valency > 2j, or the open vertex must have valency 
> 2j — 1. The minimal absolute Euler characteristic |x(r')| in the first case 
is 1 — j. Since the Euler characteristic is calculated after the open vertex is 
removed, the minimal absolute Euler characteristic in the second case is — j 
(there must be at least j edges.) Hence such graphs do not have minimal 
absolute Euler characteristic. It follows that the only labelled graph (r, £) 
with — x(r') = j — 1 with Ie(T) containing f±-*\o) is given by: 

• Qi f C t/ijwith V = 1,1 = j; the unique graph has no open vertex, 
one closed vertex and j loops at the closed vertex. The only labels 
producing the desired data are those £ p which assign all endpoints of 
all edges labelled the same index p. 

• Assuming the up/down symmetry, the corresponding Feynman ampli- 
tude I e (T) has the form 2rL(h u ) j f {2j) (0) + ■■■, where again • • • refers 
to terms with < 2j — 1 derivatives. 
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To prove the last statement about the Feynman amplitude, we first ob- 
serve that /^(0)± appears linearly in Ig (T). Indeed, the only labelled 
graphs which produce this datum give the same label to all endpoints of all 
edges, corresponding to applying only derivatives in a single variable 
Using (41), an examination of (39) shows that the coefficient of f^ 2j \0) of 
I k (a So ,C+) equals 



2r 

E ( ft 

= l,p=l 



ppy 



_d_ 

OX, 



2.7 



£ + (. . . , Xk, ■ ■ ■ X2r) — L 



=i, P =i 



/f>(0). 



To complete the proof it suffices to observe that the diagonal matrix elements 
h± are constant in p and in the sign ±. 

Had we assumed that the obstacles were both up-down and right-left 
symmetric as in [Z3], then we would already have solved the inverse spectral 
problem in this step, since all odd Taylor coefficients of / vanish at x = 
and we will be able to recover all even ones. 



6.2.2 The special term: The data fg j 1] (0) 

We now consider the more difficult odd coefficients /± 2j ' -1 ^(0), which will 
require the attributes of the amplitude (39) detailed in (40). 

We again claim that the Taylor coefficients f±~* ^(0) appear first in the 
term of order /c _J,+1 , and we enumerate the labelled graphs which have a 
combinatorial structure capable of producing /± 2,7_1 ^(0) as a factor in Ie(T) 
in the k~i +1 term. In fact, only the main term I'm^c w ^ produce such a 
term. The following is proved in [Z2]: 

Lemma 6.2. We have: 

(i) There are no labelled graphs with — x' '■= — x(r') < 3 ' — 1 f or which Ig(T) 
contains the factor /^^(O). 

(ii) There are exactly two types of labelled diagrams (T, £) with x(r') — — J + l 
such that Ie(T) contains the factor. f±~* ^(0). They are given by: 

• ^2j+i' 3 '° ^ ^2,j+i with V = 2, / = j + 1): Two closed vertices, j — 1 
loops at one closed vertex, 1 loop at the second closed vertex, one edge 
between the closed vertices; no open vertex. Labels t p>q : All labels at 
the closed vertex with valency 2j — 1 must be the same index p and all 
at the closed vertex must the be same index q. Feynman amplitude: 
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• G^j+i 3 ' C G2J+1 with V = 2,1 = j + 1; Two closed vertices, with 
j — 2 /oops at one closed vertex, and with three edges between the two 
closed vertices; no open vertex. Labels i v ^ q : All labels at the closed 
vertex with valency 2j — 1 must be the same index p and all at the 
closed vertex must the be same index q; (h±Y~ 2 (h±) 3 D 2 j~ x C±D\ C± ~ 

(/i^)^- 2 (^) 3 / (2j - 1) (o)/ (3) (o)- 

This Lemma requires a careful search through the diagrams and an ex- 
plicit calculation of the amplitudes. It is the most 'unstable' step of the 
proof, in which any calculational error could ruin the result. For instance, it 
turns out that in addition to the two diagrams which contribute, there are 
three more which have a combinatorial structure capable of contributing the 
datum / 2j ' -1 (0) ; but the corresponding amplitudes turn out to vanish due 
to special attributes (40) of the Taylor expansion of the amplitude at the 
critical point. We should note other attributes in addition to (40) are needed 
to show that the Feynman amplitude for these diagrams vanishes. 

6.2.3 Other terms 

We also have: 

LEMMA 6.3. Other terms do not contribute the data / (2j) (0), / (2j_1) (0) in 
the term of order 

The proof of this is easy. It merely requires checking how the regular- 
ization compares to the stationary phase expansion. As noted above, the 
regularized integrals are lower order yet do not involve higher derivatives of 
the amplitude or phase. 

7 Recovering the domain 

In this section, we sketch the proof of the main result, Theorem 1.1. The 
proof is virtually identical to the proof that one can determine a Z 2 -symmetric 
analytic domain from the wave invariants of a Z 2 - symmetric hyperbolic 
bouncing ball orbit. In fact, the proof requires fewer assumptions, since L 1 
is a resonance invariant (being the least non-zero singular time in the wave 
trace.) We may then inspect the wave invariants (i.e. the Balian-Bloch trace 
invariants) at the iterates Y of 7. Give the expressions above for the wave 
invariants, it suffices to analyse the Hessian coefficients. 
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7.1 Poincare map and Hessian of the length functional 
at 7 

We will need to recollect certain facts about the Hessian of the length func- 
tion at critical points corresponding to the iterates y. In most respects,the 
discussion is identical to the interior case, so we refer the reader to [Zl, Z2] 
for background on the linear Poincare map P 7 and its relation to the Hessian 
of the length function. We will only consider the Z 2 -symmetric case since it 
is the one relevant to our result. 

Henceforth we assume O is convex, so that the bouncing 7 is necessarily 
hyperbolic and the eigenvalues of its Poincare map P 7 are of the form {e ±A }. 
They are related to the geometry by 

cosh a/2 = ^1 - ^j, (43) 

where R is the common radius of curvature at the endpoints aj. For details, 
see [PS][KT]. When O is not convex, 7 could be elliptic. We refer to [Z2] for 
the calculation of the wave invariants in that case. 

We will need formulae involving matrix elements of the inverse Hessian 
matrix H^ 1 = (h pq ) at the Z 2 -symmetric bouncing ball orbit. As in [Z2], the 
Hessian in x — y coordinates is the circulant matrix 

H 2r = C(2 cosh a/2, 1,0,..., 0,1), (44) 

that is, the matrix whose rows (or columns) are obtained by cycling the first 
row (or column) (see [Zl, Z2] for background). Therefore H 2r is diagonalized 
by the finite Fourier matrix F of rank 2r (see [Zl, Z2] for its definition and 
background). Exactly as in [Zl, Z2], we have 

H 2r = FMiag (2 cos a/2 + 2, . . . , 2 cos a/2 + 2 cos ^— - ^ ] F. (45) 



An important role in the inverse problem is played by sums of powers of 
elements of columns of H^ 1 . From (45), it follows that 

= F ' ( d ' ag ( 2^/2T2 ' " ' 2 coW2 + 2co S ^ ) ) F < 46 > 
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Hence, the first row [H 2l 1 ]i = (h 11 , . . . , h l2r ) (or column) of the inverse is 
given by: 

hlq ^twr (47) 



fc=0 



where w = and where p a , r ( z ) — 2cosho;/2 + z + z n 1 , in [Zl, Z2] we will 
need that, for any p, 

y h pq = — -. (48) 

^ 2 + cosha/2 v ; 

7.2 Domain recovery 

We now prove by induction on j that the Taylor coefficients / 2j,_1 (0), / 2j, (0) 
can be determined from B^rj as r varies over r = 1,2,3, ... . At first, we 
assume / 3 (0) ^ 0. 

It suffices to separately determine the two terms 

2(^) 2 {^(0) + ^^ / (3)(0) / (-)(0)}, 

r ^ > (49) 

and ^(^) 3 \fV(0)fto-V(0). 

^ q=l > 

As discussed in [Z2], the terms decouple as r varies if and only if 
F r (coso;/2) := X^li(^2r) 3 * s non-constant in r = 1, 2, 3, . . . . 
By the explicit calculation in [Z2], we have: 

q=l 

2r \ 

kl M=o ( cosh «/2 + cos ) (cosh a/2 + cos *f ) (cos a/2 + cosh ) ' 

It is obvious that the sum is strictly increasing as r varies over even integers. 

We now begin the inductive argument. From the j = term we determine 
/"(0). Indeed, (1 - L/ (2) (0) = cos(/i)a/2 and a is a wave trace invariant. 
From the j = 2 term we recover / 3 (0), / 4 (0). The induction hypothesis is then 
that the Taylor polynomial of / of degree 2j — 2 has been recovered by the 
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j — 1st stage. By the decoupling argument we can determine f 2j (0), / 2j ' -1 (0) 
as long as / 3 (0) ^ 0. 

As discussed in [Z2], the domain can still be recovered if / 3 (0) = 0, 
but / 5 (0) 7^ 0. The argument is essentially the same. Similarly if / 3 (0) = 
/ 5 (0) = 0, but / 7 (0) 7^ 0, and so on. If all odd derivatives vanish, then one 
still recovers the domain as noted above. 

This completes the sketch of the proof of Theorem (1.1). 
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